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ABSTRACT 


This report describes a series of low-speed airfoil designs based 
on modification to the NACA 64-206 and 64 1 -212 airfoils. Designs are 
based on potential flow theory. The report describes one of a series 
of airfoil modifications carried out under Contract NAS 2-8599, Appli- 
cation of Multivariable Search Techniques to Optimal Wing Design in 
Non-Linear Flowrields. Mr. Raymond Hicks of National Aeronautics and 
Space Administration's Aeronautical Dvision, Ames Research Center, served 
as contract monitor for the study. 
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THEORETICAL EFFECT OF MODIFICATIONS TO THE UPPER SURFACE OF 
TWO NACA AIRFOILS USING SMOOTH POLYNOMIAL ADDITIONAL THICKNESS 
DISTRIBUTIONS WHICH EMPHASIZE LEADING EDGE PROFILE AND 
WHICH VARY QUADRATICALLY AT THE TRAILING EDGE 

by Antony W, Merz and Donald S. Hague 
Aevophyeica Research Corporation 

SUMMARY 

An investigation has been conducted on the Lawrence Radiation Center, 
Berkeley, CDC 7600 digital computer to determine the effects of additional 
thickness distributions to the upper surface of the NACA 64-206 and 
64j-212 airfoils. The additional thickness distribution had the form of 
a continuous mathematical function which disappears at both the leading 
edge and the trailing edge. The function behaves as a polynomial of 
order at the leading edge, and a polynomial of order at the trailing 
edge. In the present study, is a constant and e ^ is varied over a 
range of practical interest. The magnitude of the additional thickness, 
y, is a second input parameter, and the effect of varying and y on 
the aerodynamic performance of the airfoil was investigated. Results 7 
were obtained at a Mach jujs&er of 0,2 with an angle-of-attack of 6 ° on 
the basic airfoils. All calculations employ the full potential flow 
equations for two dimnesional flow. The relaxation method of Jameson is 
employed for solution of the potential flow equations. 

Earlier studies of these airfoils (References 4-6) used other types 
of additional thickness distributions. In these studies it was found that 
increases in the airfoil thickness tended to increase both the lift and 
the negative (nose-down) pitching moment. In the present investigation, 
this trend has been partially reversed, apparently because the quadratic 
curvature at the trailing edge provides a significant downward pressure 
force at the trailing edge. The result is that large reductions in the 
pitching moment can occur simultaneously with substantial increases in the 
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lift coefficient. This trend is most pronounced when the shape parameter 

IT I •' 

Cj is such as to make the leading edge very blunt, compared to the unmodified 
airfoil. 

For the range of parameters examined, the lift coefficient was nearly 
insensitive to variations in the shape parameter (less than 10%), while 
the moment coefficient showed strong sensitivity to (from 50% to 200%) . 
Increasing the thickness parameter y caused monotonic increases in lift, 
which were of smaller magnitude than those encountered in earlier studies 

fi I' 

.(References 4-6) . The moment coefficient was insensitive to changes in 
thickness for small values of and strongly sensitive to those changes 
for large values of e^. 1 \ 

Additional thickness can be made to produce significant reductions 
in peak pressure coefficient. This is particularly true for the 64-206 
airfoil, which in its unmodified form has a very high pressure peak at 
the leading edge, due to the small radius of curvature at this point. 

For both airfoils, it was found that the peak pressure can be reduced 
to much smaller values while the lift coefficient is increased, by proper 
choice of the parameters and y. Increasing the thickness first reduced 
the maximum pressure (which typically occurs at or near the leading edge) 
and then increased this pressure (which then occurs near the quarter chord) . 

A large number of parameter combinations were studied and it was 
found that the minimum pressures obtainable at a given lift are well below 
the pressures obtained in earlier studies. On the other hand, the maximum 
lift coefficient obtained with the quadratic trailing edge modification 
is some what less than the lift obtained in these studies. 

It should be noted that viscous effects are neglected in the present 
analysis. At the higher lift coefficients the effect of viscosity could 
be significant. Further investigations incorporating a viscous flow 
model are therefore desirable. 

V' 

'A 
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INTRODUCTION 



The National Aeronautics and Space Administration and others are 
currently conducting a series of theoretical and experimental studies 
to define airfoil sections having improved performance from the aspects 
of lift, drag, pitching moment, or pressure distribution characteristics, 
References 1 and 2. Analytic investigations using airfoil surface repre- 
sentations based on high-order polynomials may result in impractical 
\Yi profiles; for example, very thin trailing edge thickness distributions 

or severe reflexes in the profile. The present study employs a continuous 
polynomial arc having two free parameters, whose characteristics are 
selected to avoid such problems. Previous optimization studies using 
multivaraible search techniques. References 1, 4, 5 and 6, generally 
indicate that shape changes which provide increased lift produce unfavorable 
changes in moment characteristics. ' -Conversely, profile changes which improve 
the moment characteristics decrease the lift coefficient. 

These characteristic trends were not followed by the results of the 
present study, however, because of the strong influence of the trailing 
edge pressure on the pitching moment. On the contrary, it was found to 
be possible to reduce the magnitude of the moment and the peak pressure 
while increasing the lift of the airfoil. /. systematic variation of the 
two free parameters was carried out in order to determine the quantitative 
relationship which holds between lift, moment and peak pressure. 
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MATHEMATICAL MODELS 


Potential Flow Equation 

Potential flow analysis is based on solution of the two-dimensional 
potential flow equation 


(a 2 - u 2 ) <}. xx +'(a 2 - v 2 ) ^ - 2uv <f» xy = 0 

where is the velocity potential, u and v are the velocity components 

u = <j» x , v = 

and a is the local speed of sound determined from the energy equation and 
the stagnation speed of sound 


f a 2 = a 2 - (u 2 + v 2 ) 

Solutions are obtained by Jameson's finite difference scheme, Reference 6. 


AIRFOIL PROFILE REPRESENTATION 

Basic Airfoil 

Ordinates for the basic NACA 64-206 and 64^-212 airfoils were 
approximated by four cubic chain polynomials in the manner of Hicks 

2 3 

y. = a F. + a x + a x + a x ; j = 1,2, 3, 4 
I o. J 1. 2. 3. 

Coefficients in the four polynomial arcs are selected on the following 
basis: • H 

3=1 = Arc represents forward portion of-npper surface 

F i “ V x 

3=2 - Arc represents aft portion of upper surface 

if' p 2 = 1 

>} - ;• 

j = 3 - Arc represents forward portion of lower surface 

f 3 ■ V* 

j = 4 - Arc represents aft portion 



The coefficients a. are determined by introducing four boundary conditions 

on the airfoil profile in each of the four airfoil arcs. Crout’s method 
for triangularization and back substitution is used to solve the resulting 
system of linear equations. Note that if four points are specified on the 
aft portion (i = 2 or 4) , a discontinuity in slope occurs whore -the poly- 
nomials join. This produces a small ripple in the pressure distribution 
at thi>, juncture point. However, since the juncture occurs at a region 
of small slope (x = .5) the effect is not significant. 

Computer-generated plots of the NACA 64-206 and 64^212 airfoils, 
together with the associated pressure distributions predicted by potential 
flow theory, are shown in Figure 1. 
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Additional Thickness 


The upper surface of the basic airfoil is modified by addition of 
the thickness -distribution function, 

c i 

Ay (x) * Ax (1 - x) 


where = 2, for the present study. It is shown in Appendix A that the 
magnitude parameter. A, can be expressed in terms of the maximum thickness, 
y, by the equation, 


A = y 


2 + e. 


(2 + c.) 


Representative functions Ay(x) are shown in Figure 2, for = .25, .75 
and 1.10. The slope of the additional thickness distribution at the leading 
edge (x = 0) is infinite for < 1, and is zero for > 1, while, the 
slope at the trailing edge is zero. 


The point at which the additional thickness distribution achieves a 


maximum is given bv: . 

E 1 e 2 

measured from the airfoil leading edge. It may be noted that as e 2 increases 
the point of maximum additional thickness moves forward. 
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OPTIMIZATION STUDS US 


In previous airfoil optimisation studies, (references 4 and 5), the 
modifications to the upper surface took the form of a pair of quadratic 
arcs, which were cotangentinl at tho point x, y. These parameters arc 
respectively the chordwisc location of the maximum addition thickness 
and the value of this thickness, doth lift coefficient and moment 
coefficient were considered as performance indices in this development. 
The present study is also concerned with a two variable optimization 
problem using the leading edge thickness distribution exponent, , and 
the magnitude of additional thickness, y. 

Lift Coefficient Maximization 


In general, maximization of lift coefficient has the form 


.'here 


= Max 


4 


_ C L_ 


= fAp(x)dx 


and the integration is around the airfoil contour. The airfoil contour 
in the present study and those of references j* and 5 are completely described 
in terms of two parameters, and For the present airfoils 


* - Max [cj " Max j^C L (e 1 ,y)j B Max ^(a^ot^J 


where 


“l ~ V - “l 

1 H 

“ 2 . 1 “2 * a 2 .. 

L It 


This two variable optimization problem can be solved by use of multi- 
variable search techniques, for example, a combination of directed random- 
ray and pattern searches, references 3 and 7. 


examples illustrating this type of search procedure have previously 
been presented in references 4 and 5, However, the low dimensionality of 
the present problem (two parameters) permits the solution of optimization 
problems by inspection of graphical results. This procedure is employed 
in the present report. Other optimization problems of interest are described 
below. 
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Moment Coefficient Minimization 

Minimization of the moment coefficient has form 

* * Min [ c M ] s Min [SiM)] ^ Min [c M Co 1 *« 2 )j 

where 

C M * ^(x - 1/4) Ap(x)dx 

In previous studios moment minimization resulted in a solution directly 
opposed to lift maximization. The position of maximum thickness moved 
forward and the amount of additional thickness was minimized. Thus in those 
studies the basic airfoil had less adverse moment than ‘any airfoil generated 
by addition of the specified thickness to the upper surface of the airfoil. 

Other Optimization Criteria 

Other airfoil performance criteria can be considered, which typically 
involve compromises between lift, moment and peak pressure coefficients. 

Such modified criteria can take any of the following forms: 

1. Maximize a linear combination of lift and moment: 



2. Minimize the moment at a specific value of lift: 



3. Minimize the peak pressure at a specific value of lift: 



In the studies reported in references 4 and S, the parameters avail- 
able for airfoil modification (x and y) , both were varied over a large 
range. This permitted a straightforward interpretation of results, such 
that optimal parameter pairs for a given performance criterion could be 
determined by inspection. As noted above, this procedure is also followed 
in the present study. Free variables for the present study arc the 
parameters e 1 and y. \ 
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SYSTEMATIC AIRFOIL SHAPING 


The present study is primarily concerned with a limited but syst$t&itic 
investigation on the effects of varying tho leading edge thickness magnitude 
and distribution for the NACA 64-206 and 64^212 airfoils. The parameters 
y and are varied over ranges which are sufficient to permit qualitative 
conclusions as t0 their effects on lift, moment, and peak pressure coeffi- 
cients. For each pair of such parameters, a plot of the airfoil and its 
associated calculated pressure distribution are given in Figures 3(a) to 
. 3(o) and 4 (a) to 4 (o) . 

The first 15 of these plots relate to modifications of the 64-206 
airfoil, and the last 15 are related to modifications of the 64^-212 airfoil . 
These pressure signatures differ from those of the basic airfoils (Figure 1) 
chiefly in the magnitude of the peak pressure at the leading edge. Increases 
in hjt'ii e^ and in y tend to soften the pressure variations over the ispper 
iiitfrace, by reducing the leading edge peak and by increasing pressures over 
the central and trailing edge regions. In the case of the modified 64-206 
airfoil high values of ultimately reverse this trend and the strong 
overpressure peak reappers. 

Lift and moment coefficient variation for the two airfoils are shown 
in Figure 5, and the effects of varying the parameter and y are apparent. 
For both airfoils, the lift increases only slightly with the exponent e^, 
while the thickness y has a more pronounced influence on the lift increment. 
The adverse moment coefficient also rises sharply with additional thickness, 
while the exponent has a more significant influence on C M< Combined 
lift vs. moment results are given in Figure 6, which shows that the least 
adverse moment is obtained at a given y, or at a given C^, with the 
smallest value of e^ This corresponds to a relatively blunt leading edge 
on the airfoil. 

Variations of the peak pressure with the parameters e^ and y are shown 
in Figure 7 for the two airfoils being studied. The pressure variation 
for both airfoils is minimized at a particular y by a specific choice of 

For the 64^-212 airfoil peak overpressure at a given lift coefficient 

/j 
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is minimized by using the largest leading edge exponent value, Cj. The 

C.-C variation is more complex for the 64-206 airfoil, however, in 
Pmax 

that the constant loci cross each other. The envelope of these curves 
is given in Figure 8, and it defines the minimum peak pressure for a 
given lift. The small number of data points available in this study does 
not permit accurate cross-plotting, so the estimated minimum pressure peak 
values are shown in the cross-hatched area. Despite the uncertainties in i 
the cross-plotting procedure it is evident that the present family of 
airfoil designs produce lower peak overpressures for a given than the 
airfoils previously obtained through biquadratic modifications. 

./ ) 

Qualitative results of this parametric study are summarized in Table I . 

" a 

These conclusions follow directly from the results shown in Figures 4 to 7. 
Figure 7 also presents the minimum peak overpressures obtained with bi- 
quadratic airfoil modifications in References 4 and 5. 
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TABLE I . PARAMETRIC VARIATIONS 1*011 OPTIMIZING VARIOUS CRITERIA 



Criterion 

Exponent 

Thickness y 

Comment 


Max 

H 



Max 

Max 


64-206 

Min 

1 n [ 

\m\ 



Min 

4 

Max 

- 

Airfoil 

Min 

jc j 

1 ml 

] 


Min 

9 * 9 (c L ) 

% 


- _ 


£ l 






Min 

c 

p 

r max 

- 

c i = e i oy 

9 = 9 (c L ) 

| 










Max 




Max 

Max 

•1 :/ 
! \i 

Insensitive to 

64,-212 

Min 

Kl 



Min 

Min 


Airfoil 

Min 

h\ 



Min 

9 = y (c L ) 

“ 




£ L 




. 


Min 

c 

p 

*max 


Max 

9 = 9 (c L ) 

- 

■-> 




~ C L 


















CONCLUSION 


A numerical study has been completed of a class of modifications 
to the NACA 64-206 and 64^212 airfoils. Systematic changes in the 
upper surfaces of these airfoils were studied by independent variations 
in the thickness and leading edge thickness distribution exponent. The 
Mach number and angle-of-attack were constant during the study, and the 
results are summarized as follows: 

. 1. Pressure distribution is moderately sensitive to leading edge 
profile and to additional thickness, 

2. Lift coefficient is nearly independent of the leadin* edge profile, 
but increases with additional thickness. 

3. Adverse pitching moment increases with additional thickness and 
with increases in the leading edge profile exponent. 

u 

" 1 

4. Peak pressure for a given lift coefficient can be considerably 
reduced by careful selection of the leading edge additional 
thickness distribution exponent. Somewhat lower peak pressures 

at a given lift are possible using the present airfoil modifications, 
as compared with the '’biquadratic" modifications of References 4 
and 5. 

5. For a given lift coefficient, adverse pitching moment is minimized 
by reducing the leading edge profile exponent. 
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(c) Ej = 1.10 



FIGURE 2. UPPER SURFACE MODIFICATIONS (e 2 = 2) 
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I. MODIFIED 64-206 AIRFOIL, y = .06, Ej = .25, e 2 = 2 







FIGURE 3 (m) , MODIFIED 64-206 AIRFOIL, y = .09, Ej = .75,, c 2 = 2 
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FIGURE 3(n). MODIFIED 64-206 AIRFOIL, y = .09, e 1 = .09, e 2 " 2 
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FIGURE 3 Co).'' MODIFIED 64-206 AIRFOIL,,- 
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FIGURE 4(a). MODIFIED 64j-212 AIRFOIL, y = .03, Ej ■ .25, 
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FIGURE 4(c). MODIFIED 64j-212 AIRFOIL, y = .03, = .75, e 2 = 2 




FIGURE 4(g). MODIFIED 64^212 AIRFOIL, y = .06, e x * .5, e 2 ■ 2 

f\ If 



I: 

FIGURE 4(h). MODIFIED 64^212 AIRFOIL, y = .06, $ 1 = .75, e 2 = 2 
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FIGURE 4 (k) . MODIFIED 64^212 AIRFOIL, y = .09, e x = .25, c 2 » 2 



FIGURE 4(1). MODIFIED 64^212 AIRFOIL,? = .09, e x = .50, & 2 = 2 



FIGURE 4 (o). MODIFIED 64^212 AIRFOIL, y = .09, t 1 - 1.1, e 2 





.5 Ej 1.0 .5 c 1 1.0 

(b) 64^-212 Airfoil 

FIGURE 5. AERODYNAMIC COEFFICIENT VARIATIONS WITH EXPONENT * 

AND THICKNESS y 
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(a) 64-206 Airfoil 




Biquadratic 

Modification 


(Ref. 4) 


e 2 = 1 Modification (Ref. 6) 


2 Modification 





APPENDIX A 

CHORD WISE LOCATION OP MAXIMUM THICKNESS 


The distribution function used in this study is 

G C 

Ay(x) = Ax 1 (1 - x) 2 

und the variation of this function is smooth for 0 < x £ 1 . The 
point of maximum additional thickness occurs when 


(A-l) 


r i “1 
Cx) j^x - 


(1 - x). 


-1 


= 0 


Ay (x) = Ay (j 

This shows that the chordwise location of the point of^ maximum 
thickness is 


x = 


e l + e 2 


and the value of the maximum thickness is then found in terms of 
the parameters as 


Ay = y = A 
'max ' 


C 1 e 2 
£ 1 e 2 


(Ej + c 2 ) 


e l + e 2 


For the case studied in this report, e 0 K 1 , and the parameter 
A is 


(A-2) 


(A- 3) 


(A-4) 


A = y 


1 + e. 


(1 + C D 


(A-5) 
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